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1 Introduction 

The search for rare decays of charm particles is reaching a point where long range 
effects should be observable in the next few years. In order to perform the analysis 
of such small data sets, several techniques have become de facto standards in the last 
decade. In this talk, we will show new rare charm decay results from the FOCUS 
experiment. The analysis has been performed utilizing a bootstrap technique which 
we will describe below. 

2 Technique 

The FOCUS analysis technique has emphasized a careful approach to the treatment 
of backgrounds in a limited statistics analysis. The usual approach is to select cuts to 
optimize signal efficiency relative to background sidebands. This "blind" technique, 
where the signal region is masked off, can still lead to a downward fluctuation of the 
sidebands relative to the masked off "signal" region and a more conservative limit on 
average. 

Further, authors frequently use the technique outlined in reference [[I] to calculate 
the confidence levels used in the calculation of their limits. The approach in |]J does 
not explicitly include fluctuations in the background. Indeed, the PDG [§] suggests 
presenting a measure of the experimental sensitivity in addition to the reported limit 
whenever experiments quote a result since none of the methods suggested , includ- 
ing 0, properly deal with fluctuations in the background and bias in selecting the 
data. 

For our analysis, we chose a method Q which includes the background fluctuations 
directly into the calculation of the likelihood. The probability of finding a signal rate 
\x and a background rate b given x events in a signal region and y events in background 
sidebands is described by: 



(/i + bfe-^ (Tb) x e-W 
x\ y\ 

Where r is the ratio of the number of background events in the sideband regions 
to the signal region. The r is determined via Monte Carlo. Rolke and Lopez Q have 
shown that this method for determining an upper limit provides better coverage than 

I- ' " 

Methods presented by Rolke and Lopez have also been shown to reduce bias in 
the selection of optimal cuts 0. Briefly, one uses the computed sensitivity on an 
ensemble of bootstrapped (sample with replacement) events from the data. The cuts 
chosen by each selected set in this manner are applied to a second, independent, 
bootstrapped sample. A sensitivity as well as an upper limit are computed from the 
second sample. The sensitivity and branching ratio quoted are the median values of 
the ensemble of results from the second bootstrap. 
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Figure 1: Modes used to normalize the rare decay modes. The yields are shown for 
the loosest cut grid combination. 



The data for the analyses presented in this talk were collected using the Wide- 
band photoproduction experiment FOCUS during the 1996-1997 fixed-target run at 
Fermilab. The FOCUS detector is a large aperture, fixed-target spectrometer with 
excellent vertexing and particle identification used to measure the interactions of high 
energy photons on a segmented BeO target. The FOCUS beamline || and detector 
f?|, U have been described elsewhere. 

To provide a clean sample of D + 's and D^s, we look for D's through the 3-body 
decay chain D + (D+) —> h^/J, where the h represents a pion or a kaon. We use a 
cut grid based on kinematic variables and particle ID algorithm results that have been 
shown to be effective for other charm decays. The grid includes cuts on the significance 
of separation between the interaction and decay vertices (L/a), the confidence level 
of the vertex fit of the decay vertex, the confidence level that a particle is identified 
as a muon, the momentum of muon candidates, and Cerenkov likelihoods based on 
different particle hypotheses used to separate pions and kaons f|. The normalizing 
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modes used to compute the branching ratios for D + and are shown in Figure |l| 
for the loosest cuts in the grid. 

3 Results and Systematic Checks 

The dominant systematic effects in this analysis have been estimated to occur from 
the absolute branching ratios || used to calibrate the normalization modes. The 
systematic errors, which averaged about 7.6% for the D + modes and 27.6% for the 

modes, were added to the result using the technique outlined in @. As a check on 
the dual bootstrap, another technique was used that selected a unique cut set based 
on the results of the cut grid. The cuts used to determine the best sensitivities in the 
first bootstrap are examined for all modes and a best set is determined based on the 
most likely cut combination. This cut set is then applied to all modes once in the 
spirit of a more traditional "blind" analysis. 

The results of the analysis are presented in Table [I] below. There is good agreement 
between the dual bootstrap, the sensitivity and the single cut systematic check. Note 
that the only mode where the result and the sensitivity show a marked difference is 
in D + —> K~fi + [i + which might indicate more contamination from D + — > K~tt + tt + 
than anticipated. 



Decay Mode 


Result 


Sensitivity 


Single Cut 


D+ -» K + fi + fi- 


9.2 x 10" 6 


7.5 x 10" 6 


11.8 x 1(T 6 


D + -» K-/2+/2+ 


13.1 x 10~ 6 


4.8 x lO" 6 


12 x 10" 6 


D + -> tt+^V" 


8.8 x 10" e 


7.6 x l(T e 


7.5 x 10" 6 


D + -> 7T-fl + fl+ 


4.8 x 10" 6 


5.6 x lO" 6 


5.2 x 10~ 6 


D+ -> 


3.6 x 10" 5 


3.3 x 10" 5 


3.8 x 10" 5 




1.3 x 10" 5 


2.1 x 1(T 5 


2.0 x 10" 5 


D+ -> 7C + fl + fi- 


2.6 x 10" 5 


3.1 x 10" 5 


1.8 x 10" 5 


Df 7T>+/1 + 


2.9 x 10" 5 


2.3 x 1(T 5 


2.2 x 10" 5 



Table 1: FOCUS results with incorporated systematic errors for the modes shown. 
Each number represents a 90% upper limit for the brancing ratio of the decay mode 
listed. 



4 Comparisons to Experiment and Theory 

In Figure || below, we show the FOCUS branching ratio limits for the dual bootstrap 
technique described in the text. Our results are a substantial improvement over 
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Figure 2: FOCUS results compared to other experiments and theory. The previous 
limits, except for the E687 D + — > K~fi + fi + [12] are from Fermilab experiment E791 



|nj. The theory estimates come from JT3| (R-Parity MSSM and SM-1) and |Lj (SM- 
2). Note that the SM estimates from |H| use a formalism close to and at present 
there is some discrepancy in the invariant Mu mass behavior for the SM estimates in 
13| and \TM. The results for four-body decays measured by E791 |TTJ are not shown. 



previous results |T2|, |10j and FOCUS sets a new limit for the MSSM R-Parity violating 
prediction |]| for the branching ratio D + — > n + /j + of 8.8 x 1CT 6 @ 90% C.L.. 



I am grateful to Daniel Engh for supplying all the FOCUS results presented in 
this note. I am also grateful to Gustavo Burdman and Paul Singer for their patience 
during several very useful conversations. 
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